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ABSTRACT 

We  model  multiwavelength  afterglow  data  from  the  short  gamma-ray  burst  (GRB)  090510  using  a  combined 
leptonic-hadronic  model  of  synchrotron  radiation  from  an  adiabatic  blast  wave.  High-energy,  >100  MeV,  emission 
in  our  model  is  dominated  by  proton-synchrotron  radiation,  while  electron-synchrotron  radiation  dominates  in  the 
X-ray  and  ultraviolet  wavelengths.  The  collimation-corrected  GRB  energy,  depending  on  the  jet-break  time,  in 
this  model  could  be  as  low  as  3  x  1051  erg  but  two  orders  of  magnitude  larger  than  the  absolute  /-ray  energy. 

We  also  calculated  the  opacities  for  electron-positron  pair  production  by  /-rays  and  found  that  TeV  /-rays  from 
proton-synchrotron  radiation  can  escape  the  blast  wave  at  early  time,  and  their  detection  can  provide  evidence  of  a 
hadronic  emission  component  dominating  at  high  energies. 

Key  words:  gamma-ray  burst:  individual  (GRB  090510)  -  relativistic  processes  -  shock  waves 
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1.  INTRODUCTION 

Gamma-ray  burst  (GRB)  science  has  entered  a  new  era  with 
launch  of  the  Fermi  Gamma-ray  Space  Telescope.  The  main 
instrument,  the  Large  Area  Telescope  (LAT),  is  more  sensitive 
than  any  previous  instrument  in  the  20  MeV-300  GeV  range 
(Atwood  et  al.  2009),  whereas  the  Fermi  Gamma-ray  Burst 
Monitor  (GBM)  observes  the  whole  unocculted  sky  in  the 
8  keV^-0  MeV  range  (Meegan  et  al.  2009).  With  the  advent 
of  the  Burst  Alert  Telescope  (BAT),  X-Ray  Telescope  (XRT), 
and  UV-Optical  Telescope  (UVOT)  on  board  the  Swift  satellite 
(Gehrels  et  al.  2004),  it  is  now  possible  to  obtain  simultaneous 
multiwavelength  data  in  the  optical  to  multi-GeV  /-ray  energy 
range  from  GRBs. 

GRB  090510  is  the  first  GRB  to  provide  data  from  simul¬ 
taneous  observations  by  Fermi  (Abdo  et  al.  2010)  and  Swift 
(Hoversten  et  al.  2009)  as  well  as  by  a  couple  of  other  satel¬ 
lites.  At  a  redshift  z  =  0.903  ±  0.003  (Rau  et  al.  2009),  the 
isotropic-equivalent  /-ray  energy  release  from  this  short  GRB 
(Tgo  <  2  s)  is  Eyy so  =  (1.08  ±  0.06)  x  1053  erg  with  a  flu- 
ence  of  (5.03  ±  0.25)  x  10~5  erg  cm-2  in  the  10  keV-30  GeV 
range  (Abdo  et  al.  2010).  Fermi-LAI  detected  long-lived  emis¬ 
sion  up  to  ~200  s  after  trigger  (T0  =  00:23:00  UT,  2009  May 
10)  in  the  >100  MeV  range.  While  such  high-energy  emission, 
that  is  temporally  extended  beyond  the  keV-MeV  emission, 
was  first  detected  in  GRB  940217  by  the  Compton  Gamma-Ray 
Observatory  (Hurley  et  al.  1994),  this  feature  is  common  to  most 
GRBs  detected  with  Fermi-LKY.  Swift  XRT  and  UVOT  collected 
data  from  GRB  090510  between  ~7o  +  97  s  and  To  +  1.9  ks  be¬ 
fore  an  Earth  Occultation  (EO),  and  again  after  To  +  5.1  ks 
(Hoversten  et  al.  2009).  Swift  BAT  collected  most  data  within 
Tgo  =  0.3  ±0.1  s  (15-350  keV),  and  sparsely  between  To  +  0.4  s 
and  ~7o  +  100  s  (Hoversten  et  al.  2009). 

Smooth  temporal  evolution  of  the  flux,  F  oc  t~a,  of  the 
long-lived  emission  in  Fermi- LAT  (ay  =  1.38  ±  0.07),  and 
Swift  XRT  (ax.\  =  0.74  ±  0.03  before  EO)  and  UVOT 
(a0  j  =  — 0. SO^g1^  before  EO)  observations  strongly  suggests 
an  afterglow  origin  (De  Pasquale  et  al.  2010).  Synchrotron 
radiation  by  shock-accelerated  electrons  in  a  decelerating  GRB 
blast  wave  (Meszaros  &  Rees  1997;  Sari  et  al.  1998)  have 
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successfully  explained  much  of  the  broadband  afterglow  data 
at  radio,  optical,  and  X-ray  frequencies  in  the  pr e-Fermi  era. 
However,  fitting  combined  Fermi  and  Swift  data  from  GRB 
090510  with  simple  e-synchrotron  model  results  in  unusual 
parameter  values,  and  most  importantly  it  is  difficult  to  reconcile 
the  Fv  oc  t~av~P  temporal  relations  (De  Pasquale  et  al.  2010). 
More  complex  scenarios  have  been  proposed  to  model  GRB 
090510  data  such  as  a  radiative  fireball  in  an  e£  pair  dominated 
environment  (Ghirlanda  et  al.  2010;  Ghisellini  et  al.  2010),  an 
adiabatic  fireball  in  a  low-density  medium  and  small  magnetic 
field  (Kumar  &  Barniol  Duran  2009a,  2009b;  Gao  et  al.  2009), 
and  a  two-component  jet  (Corsi  et  al.  2010). 

Here,  we  present  a  combined  leptonic-  and  hadronic- 
afterglow  model  to  fit  multi  wavelength  data  from  GRB  090510. 
Inclusion  of  ion  acceleration  and  radiation  in  the  GRB  blast  wave 
is  a  natural  and  simple  extension  of  the  e-synchrotron  blast  wave 
model  and  has  been  discussed  by  a  number  of  authors  (Bottcher 
&  Dermer  1998;  Totani  1998;  Zhang  &  Meszaros  2001;  Wang 
et  al.  2009;  Razzaque  et  al.  2010).  We  show  that  LAT  emission 
in  the  >100  MeV  range  is  dominated  by  synchrotron  radiation 
from  protons  accelerated  in  the  external  forward  shock  of  a 
decelerating  blast  wave.  (Note  that  Razzaque  et  al.  2010  consid¬ 
ered  proton-synchrotron  radiation  from  a  coasting  blast  wave.) 
The  XRT  and  UVOT  light  curves  can  be  reasonably  reproduced 
by  synchrotron  radiation  from  electrons  accelerated  in  the  same 
external  forward  shock.  We  present  the  e-  and  ion-synchrotron 
afterglow  model  in  Section  2,  compare  this  model  with  GRB 
090510  afterglow  data  in  Section  3,  and  discuss  our  results  in 
Section  4. 

2.  SYNCHROTRON  AFTERGLOW  MODEL 

With  a  coasting  bulk  Lorentz  factor  of  To  =  103T3  and  an 
isotropic-equivalent  kinetic  energy  E^, iso  =  1 055  £55  erg,  the 
deceleration  timescale  for  an  adiabatic  blast  wave  in  a  medium 
of  uniform  density  n  cm  3  is  (Blandford  &  McKee  1976;  Sari 
et  al.  1998) 

rdec^  1-9(1 +  z)(£55/n)1/3r3 8/3  s.  (1) 

Later  the  bulk  Lorentz  factor  evolves  as 

T  «  763(1  +  z)v\E55/n)mtf3/s,  (2) 
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where  ts(>  t2cc)  is  measured  in  seconds.  At  the  deceleration  time 
t  =  tdec.  r  ~  ro/22 3/4.  The  radius  of  the  blast  wave,  R  =  4T2cf, 
is  given  by 


For  ions  of  atomic  mass  A  and  charge  Z  the  saturation  Lorentz 
factor  is  calculated  by  equating  the  acceleration  time  to  the 
shorter  of  the  dynamic  time  and  synchrotron  cooling  time  as 


R«1.4x  1017(1  +  z)~1,4(E55ts/n)l/4  cm.  (3) 

The  jet-break  time  at  which  T  ~  1  (Sari  et  al.  1999),  where 

6q  —  O.10_i  is  the  jet  opening  angle,  is  given  by 

(jet  «  105(1  +  zXFss/rc)173!?!7!3  s.  (4) 


The  fractions  of  energy  injected  in  a  forward  shock 
(Blandford  &  McKee  1976)  that  channel  into  electrons2  and 
into  ions  can  be  calculated  from  their  shock-accelerated  spec¬ 
tra.  We  assume  an  electron  injection  spectrum3  n'e(y •')  oc  y^k 

for  Ym,e  <  Ye  <  vLi.e-  Here’  Ym,e  =  Ve(m P. / me)T(t)  and  y'at  e 
are  the  minimum  and  saturation  Lorentz  factors,  respectively,  for 
the  electrons.  In  the  case  of  ions,  we  assume  an  injection  spec¬ 
trum  n'A(y'A)  oc  y'Akl  forT(t)  <  y'A  <  y'mA  and  n'A(y'A)  oc  yAkl 
for  Km, a  <  Y'a  <  Ksat,A-  Here’  Y,'„.a  =  'MHO  is  a  break  in 
the  spectrum  and  ys'at  A  is  the  saturation  ion  Lorentz  factor.  The 
fraction  of  shock  energy  carried  by  the  electrons  is 


HeVe 


k  ~  l  i  -  (YlnjYUe)'-2  . 

k-2l-(y^e/yUe)k-1’ 


2 


and  that  by  the  ions  is 


€a  —  ?AbA  - 


k,  -  1 


-2  ,  k 1- 
^2  +  r>A 


-A  + 

-1  +  rlA 


ki  ±  1 
k2  »  2. 


(5) 


(6) 


Here,  %e  and  t-A  are  the  number  fractions  of  electrons  and  ions 
that  are  accelerated  by  the  shock,  respectively,  with  an  equal 
pre-shock  number  density  n  =  ne  =  nA. 

A  fraction  eg  of  the  shock  energy  is  assumed  to  generate 
magnetic  field,  and  the  magnetic  field  behind  the  forward  shock 
is  given  by  Sari  et  al.  (1998) 


,  2.2  x  109(Z/A)4/2(£55n)1/4L1/4 

Ksat,A  «  - „  .  i  f  <  {d,A 


(1  +  z)1/40a 
_  1.2  x  1010(A/Z3/,2)f3/16 

Ks3t’A  *  (1  +  z)3/l6(f>A2elJ4(E55n3)lA6 

The  transition  takes  place  at 


;  t  ^  td,A-  (10) 


td.A  ^  1.4  X  1012(1  +2)<^e-1V7£5-55  s.  (11) 

Note  that  the  cooling  Lorentz  factor  for  ions  y'  A  — 
(A3  / Z4){mp/me)3y'  e  can  be  larger  than  the  saturation  Lorentz 
factor  (Equation  (10)). 

To  calculate  synchrotron  spectra  at  different  epoch  and  light 
curves  at  different  frequencies  arising  from  a  forward  shock,  it 
is  sufficient  to  calculate  different  spectral  break  frequencies  and 
flux  normalization  along  with  their  time  evolution  (e.g..  Sari 
et  al.  1998;  Chevalier  &  Li  2000;  Panaitescu  &  Kumar  2000; 
Granot  &  Sari  2002). 

The  characteristic  synchrotron  frequencies  for  the  electrons 
with  the  minimum,  saturation,  and  cooling  Lorentz  factors, 
respectively,  are  given  by 


hvm  e  «  7.7(1  +z)1/We(eBE55y/2t-3/2  GeV, 
hvsaUe  «  180(1  +zy5/&(p;1(E55/n)l'&t^&  GeV, 
hv(,e  «  0.5(1  +  z)“I/2e“3/2(£55«2r1/24“I/2  eV.  (12) 


A  transition  from  the  fast  cooling  (t>c  c  <  vme)  to  slow  cooling 
(vce  >  vme )  takes  place  at 


to,,  «  1.5  x  1010(1  +z)(eBrle)2nE55  s.  (13) 


In  both  the  fast-  and  slow-cooling  cases,  the  maximum  e- 
synchrotron  flux  is  given  by  (e.g..  Sari  et  al.  1998) 


B’  «  297(1  +  z)3/84/2(E55»3)1/8f73/S  G.  (7) 


In  the  regime  of  our  interest  gc,  the  Compton  parameter 

Y  =  [— 1  +  vT  +  4ee/eB]/2  —>  0  and  the  energy  loss  by  the 
electrons  is  dominated  by  synchrotron  radiation  (Sari  &  Esin 
2001). 

The  saturation  Lorentz  factor  for  electrons  is  calculated  by 
equating  the  acceleration  time  to  the  synchrotron  cooling  time 
in  the  B'  field  (Equation  (7))  as 


K sat, 


6.8  x  106  f3/16 

(l+z)3/164/2e3/4(£55«3)1/16' 


(8) 


Here,  </>“*  is  the  acceleration  efficiency  for  electrons.  The 
cooling  Lorentz  factor,  found  by  equating  the  synchrotron 
cooling  time  to  the  dynamic  time  tAyn  =  ?L/(1  +  z),  is  given 
by 

Yc,e  *  11-5(1  +zTll%e-\Elsn5Tl/itl/\  (9) 


2  Referring  to  both  electrons  and  positrons. 

3  Comoving  frame  variables  are  denoted  with  primes. 


«  52(1  +  z)~li;e(€Bn)l/2d^2E55  Jy.  (14) 

Here,  d2&/(  1028  cm)  is  the  luminosity  distance.  Note  that  the 
synchrotron  self-absorption  frequency  is  in  the  radio  band  (e.g., 
Panaitescu  &  Kumar  2000)  and  we  ignore  that  while  modeling 
optical  to  >  GeV  data. 

The  synchrotron  break  frequencies  for  the  ions  of  minimum 
and  cooling  Lorentz  factors  can  be  expressed  as  scaling  relations 
to  the  corresponding  break  frequencies  for  electrons  as 

Lh.A  —  Yj(jj a/ k\e)~ {ill e / TYl p)  Vw  e, 

vc,a  =  {A6/Z1){mp/me)5vc<e,  (15) 


and  for  the  ions  of  saturation  Lorentz  factor  (Equation  (10))  as 


^^sat  ,A 


10g3/2(»E55)3/4 
(1  +  z)i/4(p2Ats/4 


TeV;  t  <  td,A 


^sat,A  —  (-^/Z)  (jYl p I YYl e I A)V sat,e?  t  ^  td,A •  (16) 


Note  that  the  ion-synchrotron  spectrum  is  always  in  the  slow- 
cooling  regime  {vc  a  >  v,„,a)  as  opposed  to  the  e-synchrotron 
spectrum  which  can  be  in  the  fast-cooling  regime  early  and 
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Figure  1.  Modeling  of  GRB  090510  light  curves  with  proton- synchrotron  (solid 
lines)  and  electron- synchrotron  (dashed  lines)  radiation  from  an  adiabatic  blast 
wave  decelerating  in  a  uniform  density  medium.  The  model  light  curves  are 
computed  at  100  MeV  (black),  15  keV  (green),  1  keV  (red),  and  3  eV  (magenta) 
to  be  compared  with  the  LAT,  BAT,  XRT,  and  UVOT  data  points,  respectively. 
The  onset  of  the  model  light  curves  is  at  <7o  +  0.3  s  for  To  >  2400  for  a 
surrounding  medium  density  of  n  =  3  cm-3 .  The  other  model  parameters  are 
Z?k,iso  =  2  x  1055  erg,  ee  «  10-4,  €p  «  0.5,  and  e#  =  0.3.  See  the  main  text 
for  more  details. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


changes  to  the  slow-cooling  regime  later.  The  maximum  ion- 
synchrotron  flux  is 


t? max 

rv,A  ~ 


l 


$aZ3 


|Ef  +  ‘  -  1  &  A2  mi 


(17) 


for  k\  ^  1  and  k2  >  2. 


3.  MODELING  GRB  090510  AFTERGLOW  DATA 

Figure  1  shows  light  curves  at  different  energies  from  the 
combined  leptonic-hadronic  model  of  a  decelerating  adiabatic 
blast  wave  in  constant  density  medium.  With  an  initial  To  > 
2400,  the  blast  wave  decelerates  at  <0.3  s  (Equation  (1))  for  the 
parameters  used  here:  £55  2,  n  &  3  cm-3. 

We  model  LAT  emission  in  the  ~  1 00  MeV-4  GeV  range 
in  Figure  1  from  p-synchrotron  radiation  (A  =  1)  for  which 
vm  p  <  vp  <  vc  p.  In  order  to  reproduce  the  flux  decay  index 
aY  =  1.38  ±0.07  in  this  range  for  a  slow-cooling  spectrum,  one 
requires  k2  =  (4/3)ay  +  1  =  2.84  ±  0.09.  The  corresponding 
spectral  index  is  f$Y  =  (k2  —  l)/2  =  (2/3 )ay  =  0.92  ±  0.05. 
For  vp  <  vmp  <  vc  p,  the  p-synchrotron  flux  scales  as 
£„  oc  t'/2v'/3.  A  constraint  on  hvm,p  <  100  MeV  as  early 
as  &T0  +  0.4  s  requires  that  rjp  <  5  x  103  (Equation  (15)) 
with  €b  ~  0.3.  These  requirements  together  with  the  flux  level 
(Equation  (17))  needed  to  reproduce  LAT  data  constrain  the 
p-synchrotron  radiation  component.  The  fraction  of  jet  energy 
needed  in  shock-accelerated  protons  (Equation  (6))  is  ep  ~  0.5 
for  ki  <  0  and  «  10~4.  The  rise  of  the  >100  MeV  LAT  flux 
at  T  <  To  +  0.4  s  could  be  consistent  with  the  t2  rise  before 
the  blast  wave  enters  the  self-similar  regime  (Sari  1997).  The 
p-synchrotron  flux  in  the  optical  to  X-ray  is  much  below  the 
XRT  and  UVOT  data. 

We  reproduce  the  XRT  light  curve,  averaged  over  0.3-10  keV 
range,  with  decay  index  ax,  1  =  0.74  ±  0.03  before  the  EO 
at  7'  ~  7J)  +  1.43  ks  as  from  e-synchrotron  radiation.  The 
required  electron  index  is  k  =  (4/3)ax,i  +  2/3  =  1.65  ±  0.04 
for  ve  >  vm  e  >  vc  e  in  the  fast-cooling  case,  which  is  valid 


for  a  time  T  <  Ta  +  2  x  106  s  (Equation  (13)).  Note  that  the 
spectral  index,  Px,\  —  k/2  =  (2ax,i  +  l)/3  =  0.83  ±  0.02, 
is  close  to  that  of  fiy  from  p-synchrotron  radiation.  In  order  to 
produce  hvm  e  <  1  keV  at  the  beginning  of  XRT  observation  at 
T  «  To  +  100  s,  we  require  r]e  <  20 (me/mp)  (Equation  (12)). 
Together  with  parameter  w  5  x  1  0~  4  required  to  produce 
the  XRT  flux  level,  we  calculate  the  fraction  of  jet  energy  in 
electrons  to  be  €e  ~  1 0  4  (Equation  (5))  with  <pe  =  1. 

Electron-synchrotron  flux  in  the  UVOT  range  is  in  the 
frequency  range  vc,e  <ve  <  vme  and  scales  as  Fv  oc  r~<4v_1<. 
Although  the  observed  flux  fitted  with  01  o,  1  =  — 0.50/.q11i3  (De 
Pasquale  et  al.  2010)  is  different,  we  note  that  the  expected 
value  of  a  =  1/4  is  consistent  with  UVOT  data  in  the 
T  —  7o  &  600  s-1.43  ks  interval,  and  with  the  upper  limit 
atT  ^T0  + 90s. 

It  is  clear  that  the  observed  X-ray  and  UVOT  flux  decay 
indices  ax, 2  =  2.18  ±  0.1  and  a 0,2  =  1 . 1  3/(qI|I0,  respectively, 
after  the  EO  at  T  >  7|>  +  5.1  ks  (De  Pasquale  et  al.  2010)  are 
softer  than  the  e-synchrotron  emission.  If  the  jet  break  takes 
place  in  between  T  —  To  1.4  ks-5.1  ks,  then  the  expected 
decay  index  for  ve  >  vm  e  >  vc  e  is  oc  t~k  which  is  intermediate 
between  010,2  and  ax, 2  since  k  —  1 .65.  Because  of  the  idealized 
nature  of  the  afterglow  model  and  evolution  of  the  blast  wave 
during  and  after  the  jet  break  (e.g..  Sari  et  al.  1999),  the  observed 
flux  steepening  after  the  EO  could  still  be  due  to  a  jet  break. 
For  1.4  ks  <  /jet  —  To  <  5.1  ks,  the  jet  opening  angle  is 
0.16  <  9-i  <  0.26  (Equation  (4)).  If  the  jet  break  takes  place 
at  T  >  Tq  +  100  ks,  then  9- 1  >  0.8. 

As  shown  in  Figure  1,  the  BAT  flux  is  quite  noisy  and  cannot 
be  reproduced  by  either  e-  or  p-synchrotron  emission.  A  sim¬ 
ilar  conclusion  was  drawn  by  De  Pasquale  et  al.  (2010)  based 
solely  on  e-synchrotron  afterglow  model.  Sporadic  emission  in 
the  BAT  range  could  be  due  to  central  engine  activity,  working 
intermittently  at  a  much  reduced  emission  level  than  the  initial 
outburst. 

During  the  early  deceleration  phase,  the  soft  photon  density  in 
the  GRB  blast  wave  may  be  large  enough  to  induce  y  y  — >  e+e~ 
pair  production  and  photohadronic  (py )  interactions  by  protons, 
and  subsequent  cascade  formation.  The  target  photon  density 
can  be  calculated  as  n'y{e')  —  2di(l  +  z) Fv / ( R2 cT s')  from  the 
synchrotron  flux,  where  s'  =  hv’  =  hv(  1  +  z)/T. 

We  calculate  the  yy  pair  production  and  py  pion  production 
opacities  from  their  respective  cooling  timescales  and  the  dy¬ 
namic  timescale  for  the  decelerating  blast  wave  model  of  GRB 
090510.  The  opacities  for  the  y-rays  with  saturation  energies, 
both  from  the  e-  and  p-synchrotron  emission,  and  for  the  protons 
with  saturation  energies  are  shown  in  Figure  2.  The  top  panel 
shows  the  time  dependence  (from  right  to  left)  of  the  opaci¬ 
ties  at  the  saturation  energy  reached  at  that  time.  The  bottom 
panel  shows  the  opacities  versus  the  saturation  energies  reached 
within  the  same  time  interval.  Thus,  the  whole  time  interval  of 
the  top  panel  is  squeezed  to  fit  into  each  of  the  curves  in  the 
bottom  panel. 

The  saturation  energies  for  the  e-synchrotron  y-rays  scale 
with  time  as  /rvsat  e  ~  1  GeV  (Equation  (12))  for  the 

same  model  parameters  used  in  Figure  1 .  For  the  p-synchrotron, 
the  saturation  y- ray  energy  is  /jvsat,p  «  4.2 TeV 

(Equation  (16))  for  T  —  To  <  ~  3.7  x  10 l34>p  s 

(Equation  (11)).  Thus,  the  y-ray  saturation  energies  decrease 
with  time  while  the  opacities  increase  due  to  a  flux  increase  of 
the  target  photons.  However,  the  opacities  are  small  to  initiate  a 
substantial  e+e  pair  cascade  and  accompanying  radiation.  The 
same  is  true  for  photopion  cascade  (see,  however,  Asano  et  al. 
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Figure  2.  Opacities  for  the  electron-positron  pair  production  at  the  saturation 
energies  of  the  electron- synchrotron  (dashed  lines)  and  proton- synchrotron  (dot- 
dashed  lines)  emission.  Also  shown  is  the  opacity  for  photopion  ( py )  interaction 
(solid  lines).  The  horizontal  axis  in  the  top  panel  corresponds  to  the  opacities 
at  different  time,  and  in  the  bottom  panel  it  corresponds  to  the  saturation 
energies  reached  within  the  same  time  interval  for  each  lines.  We  used  the 
same  parameters  as  used  in  Figure  1  to  calculate  the  opacities. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

2009),  although  a  small  fraction  of  protons  above  Ep  >  300 
EeV  (cpp  —  1)  can  escape  as  cosmic  rays  by  converting  to  neu¬ 
trons.  It  is  interesting  to  note  that  the  saturation  proton  energy 
decreases  with  time  as  EsattP  lA\<ppxts  1/,S  EeV  because  of 

_ -3  /O 

a  decreasing  bulk  Lorentz  factor  T  ~  923 f,  '  (Equation  (2)), 
even  though  the  saturation  Lorentz  factor  increases  with  time  as 
Ksat.p  oc  t'14  (Equation  (10)). 

4.  DISCUSSION  AND  CONCLUSIONS 

We  have  fitted  the  LAT  and  XRT  light  curves  from  p-  and 
e-synchrotron  emissions,  respectively,  before  the  EO  from  an 
adiabatic  blast  wave  in  a  constant  density  medium  (Figure  1). 
The  photon  index  for  >100  MeV  emission  is  1  +  f}y  — 
1.92  ±  0.05  for  the  proton  index  ki  =  2.84  ±  0.09  and  is 
consistent  within  la  of  the  measured  value  of  the  power-law 
component  in  the  0.9  s-1.0  s  time  interval  (Abdo  et  al.  2010). 
At  earlier  time,  0.6  s-0.8  s,  where  the  power-law  component 
is  significant  and  harder,  the  agreement  is  within  3a.  However, 
emission  below  ~20  keV  in  that  time  interval,  which  we  do 
not  model,  may  contribute  to  the  hardening  of  the  power- 
law  component.  Our  model  is  compatible  with  photospheric 
emission  which  can  dominate  in  the  GBM  and  BAT  range  (e.g., 
Razzaque  &  Meszaros  2006;  Gao  et  al.  2009;  Toma  et  al.  2010). 

The  isotropic -equivalent  y-ray  energy  release  from  GRB 
090510  is  two  to  four  orders  of  magnitude  larger  than  the  typical 
range  of  1049-1051  erg  for  short  GRBs,  derived  from  pre-Fermi 
era  data  (Nakar  2007).  Our  model  requires  even  larger,  1055  erg, 
isotropic-equivalent  energy  release  in  the  jet  (see,  however, 
Freedman  &  Waxman  2001;  Eichler  &  Jontof-Hutter  2005). 
Thus,  the  jet  of  GRB  090510  must  be  strongly  beamed  for  our 
model  to  be  viable.  Depending  on  the  jet  opening  angle  derived 
in  Section  3  from  our  model  fit,  the  absolute  jet  energy  is  in 
between  ss  (3-7)  x  1051  erg  if  the  jet  break  takes  place  between 
1.4  and  5.1  ks.  For  a  jet  break  at  >100  ks,  the  absolute  energy 
release  is  >6  x  1052  erg.  Note  that  the  inferred  absolute  energy 
of  the  Fermi-LAI  GRB  sample  is  generally  high,  reaching 
1053  erg  in  some  cases  (Cenko  et  al.  2010).  We  also  note  that 
the  microphysical  parameter  ee  &  10~4  that  we  derive  for  GRB 


0905 10  using  the  combined p-  and  e-synchrotron  model,  is  much 
lower  than  the  typically  assumed  value  for  short  GRBs  in  the 
range  10" 1  to  1  0  2  which  is  derived  from  long  GRB  data  in  the 
pr e-Fermi  era  using  e-synchrotron  model  only  (Nakar  2007). 
Finally,  a  confirm  detection  of  large-amplitude  and  short-time 
variability,  synchronized  to  the  keV-MeV  emission,  in  LAT 
data  will  argue  against  their  afterglow  origin  (see,  e.g.,  Dermer 
&  Mitman  1999)  as  we  modeled  here. 

A  crucial  test  for  the  hadronic  emission  model  in  the  after¬ 
glows  of  GRBs  at  redshift  <0.5,  thus  avoiding  absorption  in 
the  extragalactic  background  light  (Razzaque  et  al.  2009;  Finke 
et  al.  2010),  may  come  from  ground-based  TeV  telescopes  such 
as  MAGIC,  VERITAS,  HESS,  and  HAWC.  An  extrapolation 
of  the  LAT  flux  in  the  >200  GeV  range,  where  e-synchrotron 
radiation  flux  is  negligible  and  internal  yy  opacity  is  small 
(Figure  2),  should  produce  similar  spectra  as  in  LAT  and  lower 
flux  from  p-synchrotron  radiation.  In  contrast,  a  synchrotron- 
self-Compton  model  (e.g.,  Dermer  et  al.  2000;  Zhang  & 
Meszaros  2001)  is  expected  to  produce  a  spectral  hardening 
and  flux  increase  in  the  TeV  range. 
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